Theoretical infra-red, Raman, and Optical spectra of the BsgNse cage 



Rajendra R. Zope^, Tunna Baruah^'"^, Mark R. Pederson^, and Brett I. Dunlap* 
^Department of Chemistry, George Washington University, Washington DC, 20052, USA 
^ Code 6392, Center for Computational Materials Science, 
US Naval Research Laboratory, Washington DC 20375, USA 
^Department of Electrical Engineering and Materials Science Research Center, 
Howard University, Washington DC 200059, USA and 
'^Code 6189, Theoretical Chemistry Section, US Naval Research Laboratory, Washington, DC 20375 

(Dated: February 2, 2008) 

The BaeNse fuUerene-like cage structure was proposed as candidate structure for the single-sheU 
boron-nitride cages observed in electron-beam irradiation experiment. We have performed all elec- 
tron density functional calculations, with large polarized Gaussian basis sets, on the BseNag cage. 
We show that the cage is energetically and vibrationally stable. The infra-red, Raman and optical 
spectra are calculated. The predicted spectra, in combination with experimentally measured spec- 
tra, will be useful in conclusive assignment of the proposed BseNse cage. The vertical and adiabatic 
ionization potentials as well as static dipole polarizability are also reported. 
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The discovery of carbon fuUerenesP], the synthesis of 
boron nitride (BN) nanotubesQ, and the fact that the 
BN pair is isoelectronic with a pair carbon atoms, led 
to the search for fuUerene analogues of BN. The car- 
bon fuUerenes are made close by introducing 12 pen- 
tagons in the hexagonal network of carbon atoms. The 
exact BN analogues of fullerenes are not preferred as 
the presence of pentagonal ring does not permit alter- 
nate sequence of B and N atoms. However, by Euler's 
theorem, it can be shown that fuUerene-like alternate 
BN cages can before, using six isolated ,,„a„BQ. 
The square contains four BN bonds with alternate boron 
and nitrogen atoms. Several theoretical works have 
reported the possible candidate structures for the BN 
cagesHQQIililiQQ. 

On the experimental side, there are also a few rep orts 
of synthesis of hollow BN structures^ Q El H Q- 
In a series of experiments by Oku and coorkwers, boron 
nitride clusters of different sizes were reported |lCllllLll2j . 
Perhaps the most interesting result of these works was 
the production of B24N24 in abundance. We performed 
theoretical calculations [siQ on B24N24 , and think that 



the roundness of the octahedral B24N24 cluster is the 
most likely explanation of its unexpected abundance. 

In another set of experiments, Stephan et al. irra- 
diated BN-powder and BN materialjl^ and observed 
small BN cage-like molecules. The irradiated deriva- 
tives were either close-packed cages or nested cages. The 
most observed cages in the experiment had diameters in 
the size range 4 to 7 A; the octahedral B12N12, BigNig, 
and B28N28 cages were proposed as possible structures. 
Subsequent tilting experiment in an electron microscope 
by Goldberg and coworkers permitted viewing of the 
irradiation-induced BN cages in different directions and 
corroborated the proposal of octahedral structures for the 
observed BN cages [l^. From the high resolution trans- 
mission electron microscopy images (HRTEM), they ob- 
served single-shelled cages of size 9-10 A, with rectangle- 
like outlines. The BagNag cage has roughly similar di- 
mension and was proposed as candidate structure by 
Alexandre e, .L ^ and OKn e, 0, Alexander 
et al. also performed a finite-range pseudopotential den- 
sity functional calculation on the BsgNsg and showed that 
the B36N36 cage is energetically stable and has large en- 
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TABLE I: The optimized coordinates of ineqivalent atoms in 
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Y 
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-1.340 


4.713 


4.713 
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-7.771 


-1.224 


-1.224 
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6.099 


2.739 


2.739 
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1.457 


4.876 


4.876 


N 


8.189 


1.504 


1.504 


N 


-6.057 


-2.702 


-2.702 



ergy gap between its highest molecular orbital and low- 
est molecular orbital (HOMO-LUMO). These conclusions 
were subsequently confirmed by all electron calculations 
using the Slater- Roothan (SR) method^. These calcula- 
tions show that the proposed BaeNag cage has overall size 
of about 8 A, smaller than single-shelled BN cages ob- 
served in the experiment. So far the proposal of BsgNag 
as the observed structure has not yet been verified in al- 
ternative experiments. It is the purpose of the present 
work to ascertain the stability of proposed cage structure 
of B36N36 with respect to small distortions and to pro- 
vide its spectroscopic data : the infra-red (IR) spectra, 
Raman and optical spectra, for comparison with exper- 
imental spectra. The comparison of the calculated and 
experimentally measured spectra has played an impor- 
tant role in structure determination in past. We hope 
that the present work will stimulate further experimental 
works, which combined with the predicted spectra, would 
unambiguously confirm the proposal of B36N36 fuUerene- 
like cage. 

Our calculations were carried out at the all electron 
level using the NRLMOL suite of codes The 
exchange-correlation effects were treated at the level 
of the generalized gradient approximation using the 
Perdew-Burke-Ernzerhof sche me( PBE-GGA) 18|. Large 
polarized Gaussian basis sets[l9j optimized for density 
functional calculations (PBE-GGA) are employed. The 
B-atom basis set consists of 5 s— type, 4 p— type, and 
3 d— type orbitals each of which were constructed from 
a linear combination of 12 primitive Gaussians. For N 
atoms, 5 s— type, 4 p— type, and 3 d— type functions com- 
prised of 13 Gaussian are used. The geometry optimiza- 




FIG. 1: (Color online) Two different views of optimized 
B36N36 cages. In the second view the BaeNae cage is ro- 
tated by 45° around a vertical axis passing through center of 
square. 

tion was performed using the limited-memory Broydcn- 
Fletcher-Goldfarb-Shanno scheme until the forces on 
each atom were less than 0.001 a.u. The self-consistent 
calculation was iterated until the energy difference of suc- 
cessive iteration was smaller 0.0000001 a.u. The vibra- 
tional frequencies were calculated by diagonalizing a dy- 
namical matrix constructed by displacing the atoms by 
small amount (2^. The IR and Raman frequencies were 
determined from the derivative of the dipole moment and 
the polarizability tensor. 

The optimization of B36N36 cage (See Fig. 0) was 
started using its optimized geometry obtained by the SR 
methodQ. The cage has Td symmetry and can be gen- 
erated from six inequivalent atoms, three for boron, and 
three for nitrogen, using symmetry operations. The op- 
timized coordinates of the inequivalent atoms are given 
in Table |l| These coordinates are with respect to the 
origin at the center of mass. The average spherical ra- 
dius of the cage, obtained as the mean of distance of each 

atom from the center of mass, is 3.94 A. The SR method 

n 

predicts a cage of 3.76 A radius Thus, the cage has 
roughly dimension of 8 A, somewhat smaller than the 
observed dimensions of cage structure in the experiment. 
Part of this discrepancy could be due to the thermal ef- 
fects as our calculations are at T = QK. The cage is 
energetically stable with respect to isolated atoms and 
has a binding energy (BE) of 8.43 eV per atom. The SR 
method gives a binding energy of 7.53 eV/atom using 
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FIG. 2: (Color online) The infra-red (top) and Raman (lower 
panel) spectrum of BseNse cage. 



the 6-311G* basis. Th present BE is higher than that 
of B24N24 cages, which have energies of about 8.3 eV £ |. 
This result agrees with our predictions by SR method 5 1 . 
The trend of larger BN cages being energetically more 
stable than smaller ones is similar to that observed in 
case of carbon fullerenes[2lt. The present calculations, 
consistent with earlier predictions by lower theories, show 
that the BsgNsg cage is characterized by a large HOMO- 
LUMO gap of 5.0 eV. The HOMO level has t2 symme- 
try. The vertical ionization potential (VIP) is calculated 
as the self-consistent energy differences of the cage and 
its positive ion with the same ionic configurations. The 
calculated VIP is 8.2 eV. This agrees well with 8.4 eV, 
obtained by the SR method^]. We have also determined 
the adiabatic ionization potential from the energy differ- 
ences of the neutral BagNag cage and its optimized singly 



charged cation. The adiabatic ionization potential is 8.1 
eV, which is an overestimate as the cation was optimized 
under symmetry constraint. The mean polarizability of 
B36N36 calculated from the finite field method is 78 A'^. 
The nuclear frame of BagNag was assumed to be frozen 
during the polarizability calculation. 

The vibrational analysis within the harmonic approx- 
imation indicate that all frequencies are real. Thus, the 
BagNag cage is vibrationally stable and corresponds to 
a minimum on the potential surface. The calculated IR 
spectrum is presented in Fig. [3 The IR spectrum is 
broadened by 6 cm~^ . The IR spectra shows two con- 
spicuous peaks at 1432 cm~^ and 1456 cm~^ ; the inten- 
sity of the latter peak (42 Debye^/amu/A^) is half of the 
former peak. A weak but noticeable absorption (inten- 
sity 6 Debye^/amu/A^) is found at 797 cm~^ . All the 
peaks have T2 symmetry. The Raman spectrum is also 
shown in Fig. [3 The Raman spectrum is characterized 
by multiple peaks in the frequency range 150 — 900 cm^^ 
. Of the two prominent peaks in the low frequency range, 
the first one is at 165 cm^^ and is of E type. The second 
one at 188 cm~^ is due to the triply degenerate mode 
of T2 symmetry. All other Raman active frequencies: at 
354, 439, 796, 847, and 869 cm"! are of Ai symmetry. 
The absorption at 354 cm~^ corresponds to symmetric 
stretching or breathing mode. All vibrational frequen- 
cies along with their symmetry type are listed in Table 

im 

The density of states (DOS) and optical spectra of 
the BagNag cluster calculated within DFT are shown in 
Fig. H The HOMO-LUMO gap of this cluster is 5.0 eV. 
The HOMO has a predominant N character whereas the 
LUMO is delocalized over the whole cluster. The states 
near the HOMO and the LUMO have more N character 
than B character. The HOMO has t2 symmetry and the 
LUMO has ai symmetry The HOMO-LUMO transition 
is weak though symmetry allowed. The optical spectra 
shows a sharp peak at 5.4 eV which occurs from tran- 
sitions from occupied states within 0.59 eV of HOMO 
to the unoccupied states within 0.43 eV of LUMO. This 
peak has contributions from the t2 — > ^2 and e — > ^2 tran- 
sitions. It also has weak contributions from t2 — s- e tran- 
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TABLE II: The frequencies, symmetry, IR and Raman activeness of the vibrational modes. 
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FIG. 3: The density of states of B36N36 cluster. 
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FIG. 4: The optical absorption spectra of BaeNse cluster. 

sitions. There are other higher energy peaks seen at 9.7 
eV and 13.8 eV. The small peak at 9.7 eV arises mainly 
from a t2 — » ^2 transition involving occupied states lying 
3.4 eV below HOMO and unoccupied states 1.3 eV above 
the LUMO. The broad peak at 13.8 eV is due to a large 
number of weak transitions involving a large number of 
low lying states rather than any sharp strong transition. 



It should be noted that the present calculations are per- 
formed within the PBE-GGA at zero temperature and 
hence the predicted optical spectra will be useful only 
for qualitative comparison. More accurate spectra can 
be obtained by the TDDFT or GW approximation. 

In summary, density functional theory at the level of 
PBE-GGA is used to examine the vibrational stabil- 
ity of the proposed fuUerene-like cage structure of the 
BagNsg observed in a recent electron-radiation measure- 
ment. The cage structure was proposed on the criteria 
that it satisfied the isolated six square rule for alternate 
BN cages and has dimension similar to the single-shell 
BN cages observed in the HRTEM images. The calcu- 
lations indicate somewhat smaller dimension 8 A) of 
the BagNag cage than that those observed in experiments 
(9-10 A). The cage is energetically and vibrationally sta- 
ble. It has vertical and adiatbatic ionization potential of 
8.2 and 8 eV, respectively. Its mean static dipole polar- 
izability is 79 A^. We have predicted IR, Raman and 
optical spectra at the DFT GGA level. We hope that 
the predicted spectra in this work will stimulate exper- 
imental measurements of these spectra. The measured 
spectra, in combination with the predicated spectra in 
this work, will play a decisive role in confirmation of the 
proposed candidate structure. 
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